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Reconstructing the Primordial Spectrum with CMB Temperature and Polarization 
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We develop a new method to reconstruct the power spectrum of primordial curvature pertur- 
bations, P{k), by using both the temperature and polarization spectra of the cosmic microwave 
background (CMB). We test this method using several mock primordial spectra having non-trivial 
features including the one with an oscillatory component, and find that the spectrum can be recon- 
structed with a few percent accuracy by an iterative procedure in an ideal situation in which there 
is no observational error in the CMB data. In particular, although the previous "cosmic inversion" 
method, which used only the temperature fluctuations, suffered from large numerical errors around 
some specific values of k that correspond to nodes in a transfer function, these errors are found to 
disappear almost completely in the new method. 

PACS numbers: 



I. INTRODUCTION 



The cosmic microwave background (CMB) anisotropics contain important pieces of information on physics of the 
' early universe. The recent precise data of the Wilkinson Microwave Anisotropy Probe (WMAP) tell us that our 
universe is consistent with a spatially flat universe dominated by a cosmological constant (A) and cold dark matter 
' (CDM), the so-called ACDM model, with slow-roll inflation at its early stage which has generated Gaussian, adiabatic, 
. and nearly scale- invariant primordial fluctuations P, S S 13 ■ 

' Nevertheless, it is reported that the observed CMB angular power spectrum may have some non-trivial features 
, such as lack of power on large scales [67j , running of the spectral index, and oscillatory behaviors of the spectrum on 

■ intermediate scales. To explain these non-trivial features, many authors focused on the primordial power spectrum of 
; the curvature perturbatioii^P(fc)^nd_p^ inflation models [lll[ll[ll|llli3|llEllIl[IlllllIl 

O i' IMIIlllllllEllHllllla, HI ill IM El HI m nil . However, equally important is to understand the implication 
I ' of these observed possible non-trivial features on the primordial curvature spectrum without a theoretical prejudice. 
In this respect, there have been several att emp ts to reconstruct the primordial spectrum using the WMAP data by 
model-independent methods [sl IsIIstL Isl Isa liot . 
^ In our previous w ork, we developed a method to reconstruct the primordial spectrum directly from the observed 

CMB anisotropy [4ll[4^ . which we call the cosmic inversion method, and applied it to the WMAP flrst-year data [isf . 
Compared with other reconstruction methods such as the binning, wavelet band powers, and direct wavelet expansion 
method, we have shown that our method can reproduce fine features in P{k) with a resolution of Ak ~ 3.7 x 

■ 10-''Mpc"^ which roughly corresponds to A£ ~ 5 in the angular power spectrum, Ci. Performing a statistical 
analysis of the reconstructed P{k), we have found that there are some possible deviations from scale-invariance 
around /fc ~ 1.5 x lO^^Mpc"^ and 2.6 x lO^^^pc"^ 43). 

The method we developed, however, used only the temperature-temperature (TT) angular power spectrum. The 
CMB anisotropy contains another important information, the polarization. During the recombination epoch the 
linear polarization of the CMB is generated by the Thomson scattering of the local quadrupole component of the 
temperature anisotropy 44]. From the perspective of a cross check, it is important to reconstruct P{k) taking the 
CMB polarization into account. The CMB polarization was detected by the degree angular scale interferometer 
(DASI) for the first time ji^ and more precise data of the temperature-polarization (TE) angular power spectrum 
were released by the WMAP recently |4g. In the future, much more precise observation of the polarization will be 
carried out by the Planck satellite |68l| . 

In this paper, we present an improved cosmic inversion method that takes account of the CMB polarization spec- 
trum, and test our new method for various shapes of P{k). That is, we calculate the CMB temperature and polarization 



* Electronic address: |kogo@3fukawa. kyoto-u . ac . j p | 
t Electronic address: imTsao^yuk^al{yot^u!ac! jp] 

fEIectronic address: yokoyama@vega.ess.sci.osaka-u.ac.jp| 



2 



spectra for a given P{k) and reconstruct it from thus obtained CMB data. As a first step, we neglect observational 
errors in numerical tests and assume that the CMB temperature and polarization spectra are completely known. We 
calculate the CMB power spectra based on CMBFAST but we adopt a much finer resolution than the original one 
in both k and £. Our inversion method is an iterative method based on approximate formulas for CMB temperature 
and polarization spectra and correction factors that adjust the errors in the approximate formulas. 

This paper is organized as follows. In Sec^I we review the basic theory of the CMB polarization. In Sec. IIIII we 
extend our cosmic inversion method and propose a new inversion formula that uses both the CMB temperature and 
polarization spectra. In Sec. II VI we test our new method. We find that the primordial spectrum can be reconstructed 
with much better accuracy. In Sec. we propose the idea of constraining the cosmological parameters based on our 
new method. Finally, Sec. IVIl is devoted to conclusion. 



II. BASIC THEORY 

Here we summarize basic equations for the CMB temperature fluctuations and polarization based on [4^ liM. , 
and their approximate formulas which are used in our inversion method. We deal only with scalar-type perturbations 
and assume a spatially flat universe with Gaussian and adiabatic primordial fluctuations. We leave inclusion of 
tensor-type perturbations for future work. 

It is convenient to describe polarized radiation by the Stokes parameters [Hfl l5ll |. Consider a monochromatic 
electromagnetic wave propagating in the z-dircction with an angular frequency loq. The components of the electric 
field are written as 

Ex = 8x{t) cos[LJot - ipx{t)], (1) 

Ey = £y{t) COS[LJ0t - ipy{t)]. (2) 

The Stokes parameters are defined as 

/ = {£j+£y^), (3) 

Q = {£.^-£y^), (4) 

U = {2£x£y cos{ipx - ipy)) , (5) 

V = {2£x£ysiii{ipx - ipy)), (6) 

where the brackets denote the time average over the modulation of the amplitudes and the phases. Since the Thomson 
scattering generates only linear polarization, we can neglect V which describes the degree of circular polarization. 
Under a rotation of an angle ip on the {x, ?/)-plane, / is invariant while Q and U are transformed as 

Q' ^ Qcos2il; + Usin2'ip, (7) 

U' = -Qsin2V' + C/cos2V', (8) 

which are identical to 

iQ±iUy ^e^^'^iQ±iU). (9) 

This means that {Q ± ill) have spins ±2, respectively. 

We denote the temperature fluctuations and the polarization of radiation at the conformal time t] and the comoving 
spatial position x propagating in the direction n as 0(ry,a;,n) and (Q ± iU){ri,x,n), respectively. In general, these 
are expanded using the normal modes defined as 



,G,„^Hr\/^i^,y,,„(n)e*-, (10) 



where k is the comoving wavenumber and sYim is the spin-s weighted spherical harmonic function |52l l53l l54j. For 
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each Fourier mode we can work in the coordinate system where k \\ z. In the case of scalar perturbations, m = 0, 



e(ry,a;,n) 



(27r)= 
(Pk 



d^k 



(7+2)! 



{Ei±tBi){rj,k)Pi{fi), 



(11) 



(12) 



where n = k ■ n and is the associated Legendre function, P™(/i) with m — 2. Here and iSf represent the 
so-caUed E-mode and B-mode polarizations with parities (—1)^ and (— l)^"*"^, respectively. Note that for scalar-type 
perturbations, vanishes because of its parity. Using these multipole moments, the angular power spectrum is 
expressed as 



2£+ 1 
An 



27r2 



k 2i+l 



(13) 



where X and X are either {— T) or E. The Boltzmann equations for each Oeiv^ k) and E£{r], k) can be transformed 
into the following integral form {t> 2): 



2£+l 



Eeivo,k) 
21+1 



drj 



(Oo + *)V(r;) + (* - $)e-"('')l ji{k^i^) 



VtV{v)f,{kATj) + -n2V{7j) [3j;'(fcA7?) + j,(fcA77)] }, 



{£-2)1 



(fcAr/) 



2 ' 



(14) 
(15) 



where 112 = (©2 ^ \/6£'2)/10, Ar/ = rjo — r] with r/o being the conformal time at present, and the overdot denotes a 
derivative with respect to the conformal time. Here Vb is the baryon fluid velocity, ^' and $ are the Newton potential 
and the spatial curvature perturbation in the Newton gauge, respectively ||55|, and 



EE fe-^^") 



r(7?) 



TdT], 



(16) 



are the visibility function and the optical depth for Thomson scattering, respectively. In the limit that the thickness 
of the last scattering surface (LSS) is negligible, we have Vj-n) y Sljj — rj^) and e^'^^''^ ~ Oirj — 77*), where 77, is the 
recombination time when the visibility function is maximum 56] . 

As we did in our previous papers |42l l43j| . we now take into account the thickness of the LSS approximately to 
perform the integrals in Eq. p5|) as well as in Eq. p4|l analytically. For the polarization, this is crucial since the CMB 
polarization is mainly generated within the thickness of the LSS. The approximation is to neglect the oscillations of 
the spherical Bessel functions in the integrals. Applying this approximation to Eq. H14|l and H15|> . we have 



Qf(?7o,fe) 
2^+1 



Eiivo,k) 

2e + i 



'(.end 

drj 

»?.ond 



(60 + *)V(ry) + (* - $)e-" 



iv) 



je{kd) 



+ 



dr,Q,{r,,k)V{Ti) \ j[{kd) = 



2^+1 



H^ + 2)! 
(^-2)! 



drj 



Ho 



?7*start 



(fcA77)2 



v{ii)\n{kd) 



2£+l 



(17) 
(18) 



where d = r]Q — rj^ is the conformal distance from the present to the LSS and ?7*start and ry^end are the times when the 
recombination starts and ends, respectively. We have also replaced VJ, by Oi and neglected the quadrupole term in 
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Eq. (|14() . by adopting the tight coupling approximation [H^. The transfer functions, f{k), g{k), and h{k) are defined 
by 



(eo + *)(77,A:)V(7/) + (^'-$)(r;,fc)e-^(")J = /(fc)$(0,fe), 

If.ond 



d?7ei(r7,fc)V(77) = <?(fc)$(0,fc), 

IJ.start 



(19) 
(20) 
(21) 



Given the cosmological parameters, one can calculate transfer functions numerically, e.g., by using CMBFAST. The 
results are depicted in Fig. ^ 

The primordial power spectrum of the curvature perturbation is defined as 

F(fc)^(|$(0,fe)p). (22) 

Substituting Eqs. lfT7|l and (fT^ into Eq. ifT^ . we obtain the approximated TT, EE, and TE angular power spectra. 



2^+1 ^TT, app 
2^+1 ^EE. app 

2^+1 r^TE, app 



2^+1 



27r2 

2^+1 (^ + 2)! 
27r2 (^-2)! 



^fc3p(fc) [/(fc)j,(fcd)+g(fc)j;(fcd)]^ 



^k^p{k) [h{k)n{kd)f 



-c 



2£+l {£ + 2)1 



27r2 



'dk 



i£-2y.J, k 



k^'Pik) [f{k)u{kd)+g{k)j[{kd)] h{k)Mkd). 



(23) 
(24) 

(25) 



III. INVERSION METHOD 



Here we describe our method of the reconstruction of P{k). Let us first briefly review the case of the temperature 
fluctuations proposed in 0, • The CMB temperature fluctuations can be quantifled by the angular correlation 
function defined as 

°° 2/ + 1 

C^T{0) = {e{hi)Q{n2)) —^Cj^Piicosd), cos6' = ni-n2. (26) 

^ — ' 47r 

Here we introduce a new variable r instead of 6 defined as 

Q 

r = 2dsin-, (27) 

which is the conformal distance between two points on the LSS. Substituting Eq. H23|l into Eq. H26|l and using the 
formula: 

oo . , 

Y,i'2i + i)Pi{cose)n\kd) = (28) 

we can derive formulas for the approximated angular correlation functions in terms of P{k). In the small-scale limit 
r d, its derivative is written as 

-^{r3C^^'^PP(r)} = — / dfcF(fc) [/2(fc)fc2rcosfcr + {2/2(fc)+g2(fc)}fcsinfcr] . (29) 
r or 2'K'^ Jq 

Integrating by parts and applying the Fourier sine formula, we obtain a first-order differential equation for P{k), 

-ef{k)P'{k) + [~2fc2/(fc)/'(fc) + kg\k)] P{k) 

= 4^/ dr - {r^ C^^' (r)} sin kr = S^^ik). (30) 
Jo rdr 
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This is the basic equation for the inversion of the TT angular power spectrum to the primordial curvature perturbation 
spectrum. Note that f{k) is an oscillatory function and the above differential equation is singular at f{k) — 0. 
However, this can be regarded as an advantage. Namely, we can find the values of P{k) at zero points of f{k), say 

P(^^^-ir^) for /(M = o, (31) 

assuming that P'{k) is finite at = kg- Thus, without solving the differential equation, it is possible to know a rough 
overall feature of the spectrum. We can then solve Eq. (|30|1 as a boundary value problem between the singularities. 
Now we consider the case of the polarization. Like Eq. (|26|l . we introduce the following quantities: 

^""W - E Cf-P^icose)^ (32) 



^..(,) ^ t'-^J^{^CrP.i^ose). (33) 



1=0 



These are defined in such a way that the factors (^ + 2) !/ (^ - 2) ! and ^(^ + 2)!/(^ - 2)! in Cf ^ and Cj^ , respectively, 
are canceled out so that we can use the formula (|28|l . Note that they are not the conventional angular correlation 
functions but just quantities that are convenient for inversion. Substituting Eqs. 124|l and H25|l into Eqs. H32|) and 
()33|l . respectively, in the small-scale limit r <g; d, we have 

1 1"°° 

rC^E, ^pp,^ = — / dkkh^{k)P(k)smkr, (34) 
27r^ Jo 

^(^Ti3, app/N ^ / dkkf(k)h(k)P(k)smkr. (35) 

27r^ Jo 

We obtain expressions much simpler than that for the TT case. Applying the Fourier sine formula, we obtain the 
algebraic equations for P(fc), 

/•oo 

kh'^{k)P{k) ^ An drrC^^'''PP{r)smkr= S'^'^ik), (36) 
Jo 

kf{k)h{k)P{k) ^ 4tt drrC'^^'^PP(r)sinfcr = 5'^^(fc). (37) 
Jo 

In this case, we can find P{k) at the values of k where h{k) ^ for EE, and at those where f{k) ^ and h{k) ^ 
for TE without solving a differential equation. The equations (|36ll and (|37|l are the basic inversion formulas for the 
EE and TE angular power spectra, respectively. 

In practice, however, we find that numerical errors become large around the singularities of Eq. H36(l or 1)37(1 . where 
the prefactor of P{k) vanishes, if we use either of them for inversion. We encountered a similar numerical problem in 
the TT case when we used Eq. (|3U|I for inversion, particularly when observational errors in the CMB power spectrum 
were taken into account |43j| . Here, as a remedy for this numerical problem, we propose a new method that combines 
the TT and EE formulas. Since the zero points of f{k) and h{k) are different from each other as shown in Fig.^ we 
expect this to suppress numerical errors around the zero points of both /(fc) and h(k). 

Multiplying Eq. ^ by some factor a (which we take to be independent of k for simplicity), and adding it to 
Eq. l|3Ufl . we obtain the combined inversion formula as 

- k^f{k)P'{k) + [-2k^f{k)f'{k) + kg^{k) + akh^ik)] P{k) = S'^^{k) + aS^^{k), (38) 

and the boundary conditions are similarly given by the values of P{k) at zero points of f{k) as 

^'w^Qt^^SSl ^<«^»' '''' 

assuming that P'{k) is finite at fc = kg. The contribution of EE is controlled by the parameter a. If we take an 
appropriate value of a so that the contribution of EE is comparable to that of TT, the solution of Eq. (|38|l becomes 
numerically stable even around the singularities because the contribution of EE dominates near the singularities of 
TT given by f{k) = 0, and vice versa. We investigate this in the next section. 
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We do not use the TE formula, Eq. H37|l . because it is singular not only at h{k) =0 but also at /(fc) = where the 
TT formula is also singular Q| . Perhaps the TE spectrum can be used as a consistency check when we have accurate 
temperature and polarization maps at hand. However, here we do not discuss it. 

For either TT or EE, the approximated spectrum (j^-^'^^pp (^X — T or E) used to obtain the inversion formula 
has relative errors as large as about 20 — 30% when compared with the exact one Cf'^''^^. We show Cf^'°^ and 

course, since the observed spectrum '° is to be compared with 
it is necessary to correct the errors. For this purpose, we first define the ratio, 

I^XX. cx 

bfx ^ ^ (40) 

^ /-iXX, app ^ ' 

Interestingly, this ratio is found to be almost independent of P(k) not only for TT but also EE. Figure O shows 
hj'^ and bf^ for a scale-invariant P{k) and a spectrum with a peak and a dip. In this case, the difference is found to 
be less than about 5%. We use this property of hf^ to perform an iteration of inversion as previously done in the 
TT case El El. 



XX fO) 

For self-containedness, let us recapitulate the iterative procedure. First we calculate = 
^xx, cx(o)^^xx, app(o)^ for a fiducial primordial spectrum P^^\k) such as a scale- invariant one. Then we divide 

the observed spectrum by bf^''"^^; Cf'^'°^^/bf'^''"^\ This gives a better guess to Cf^'^'^'^. We use it to 

evaluate S^^{k) in the right-hand side of Eq. and solve this equation. Let us denote the reconstructed pri- 
mordial spectrum from Cf'^'°^^/bf'^'''^^ by P''^\k). Then we can repeat the same procedure and obtain P'''^\k) 

from Cf^'°^^/bf^'^-^\ where bf^''"^^ is calculated from p(^^(/c). Repeating this procedure as many times as it is 
necessary, we can find P{k) with much better accuracy. The n-th iteration is summarized as 

^XX.cx(n-l) ^XX,ohs 

■ ^ p'^-'Hk) - bf--'--''-^^kx:,^ - ^ p^-Hk) (41) 

We confirm that this iteration converges within a few times as discussed in the next section. 

IV. NUMERICAL TESTS 

To test our inversion method including the CMB polarization, we perform the reconstruction of P(k) from Cf^^ 
for several P(k) whose shapes are given by hand. As mentioned before, we assume here an ideal situation in which 
the CMB data contain no observational error. We use Cf-^ up to ^max = 1500 for both TT and EE. As a fiducial 
spectrum, P^°\k), for the iteration, we choose a scale-invariant spectrum k^P{k) = const., and assume the cosmo- 
logical parameters as ft, = 0.70, fib = 0.050, f^A = 0.70, ri„i = 0.30, and r = 0.20. In this case, the positions of the 
singularities for f{k) = are at kd ~ 70, 430, 680, 1030, • • • , and for h{k) = are at kd ~ 230, 560, 860, 1180, 
• • • , where d ~ 1.36 x lO-^Mpc. We solve Eq. ^ as a boundary value problem between the first and fourth TT 
singularities for f{k) = 0. Hence the range of the reconstructed P{k) is 70 < kd < 1030, where kd roughly corresponds 
to £. When proceeding from the {n — l)-th to n-th iteration, we need p("-i)(fc) up to kmaxd — 2£niax to Calculate 
^xx, (n-i)^ Hence we assume the scale-invariance outside of the reconstructed range. 

First we examine the dependence of the reconstructed P{k) on the parameter a appeared in Eq. H38I) . which controls 
the contribution of EE, assuming that the cosmological parameters are known. We assume that P{k) has a peak and 
a dip expressed as 



k^P{k) = A <^ 1 + ai exp 



1 + 02 exp 



{k-k^y 



^2^ 



(42) 



and we take oi = 03 1.0, fci = 0.03Mpc"\ fca = 0.06Mpc"\ cti = 0.01Mpc~\ and (72 = 0.005 Mpc"\ FigureH 
depicts the reconstructed P{k) with four different values of a. The case a = corresponds to the result of the original 
cosmic inversion method. We show only P'^^^(fc) for each a, which means that no iteration is performed, to see the 
effect of the contribution of EE explicitly. We find that the reconstructed P{k) is in a good agreement with the 
original one even near the singularities for a = 10^^, while spurious sharp features appear near the TT singularities 
for a much smaller than 10^^ and those appear near the EE singularities for a much larger than 10^^. This behavior 
is reasonable because the contribution of EE is found to be of the same order of TT for a ~ 10"'^^"^^ by comparing 
their magnitudes in the source terms of Eq. H38|) . Therefore, we adopt a = 10^^ in the following reconstructions. It 
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may be noted that the origin of this number, lO'^'*"^^, comes dominantly from the fact that the transfer function 
h{k) for EE is intrinsically smaller than the transfer functions f{k) and g{k) for TT by a factor ~ 10"^ (see Fig. ^ 
and their squares are contained in the left-hand side of Eq. . 

We perform the inversion for several shapes of P{k), assuming that the cosmological parameters are known. 

(i) Spectrum with a peak and a dip 

Figure m shows the results of the iteration for P{k) given by Eq. H42() . We see that the iteration converges within 
a few times and the relative errors are reduced to about a few percent level. In it has been shown that 
in the case of using the TT spectrum only, which is identical to a = here, one can reconstruct P{k) with 
about a few percent errors by iteration, cutting off spurious sharp features appeared near the singularities and 
interpolating the spectrum smoothly at each step of iteration. But with our new method, it became unnecessary 
to perform this artificial modification of the spectrum during iteration. 

(ii) Spectrum with a running spectral index 

Figure shows the results for P{k) with a running spectral index given by [57| 

k \ ^ , , dn , / k 



k^Pik)=A{- , «W--(M + ^ln(^). (43) 



We take ko = 0.05Mpc"\ n{ko) = 0.8, and dn/dlnk = -0.05. 

(iii) Spectrum with an oscillatory component 

Figure shows the results for P{k) with an oscillatory component given by 



k^P{k) = A 



1 + flo sin 



(44) 



This form is motivated by possible signatures of trans-Planckian physics [581 l59l |60|. We take a = 0.1 and 

ko = 5x lO-'^Mpc"^ 

The results are quite impressive because we can reconstruct not only a global shape such as a running spectral 
index, but also a fine feature such as a small oscillation as given by Eq. (|44ll with only a few percent errors. 



V. NOVEL METHOD TO CONSTRAIN COSMOLOGICAL PARAMETERS 



We also examine how the shape of reconstructed P{k) changes if we adopt the values of the cosmological parameters 
different from the assumed values, namely, those used in making a mock sample of Cf^^ . We vary each cosmological 
parameter in the range of about 2a of its current best-fit value, that is, 0.60 < h < 0.80, 0.040 < f^b < 0.060, and 
0.60 < < 0.80, respectively, with the others fixed in the case of a scale-invariant P{k). We do not vary t since it 
affects only the normalization except for the spectrum on large scales. 

The resultant P^^) (fc) are shown in Fig. |S1 We see that if we use incorrect values of the cosmological parameters 
in our reconstruction procedure, the reconstructed P{k) is deformed from its real shape near the EE singularities for 
a — 10^^. This is because the contribution of EE is slightly larger than that of TT. If we adopt a smaller value of 
a so that the contribution of TT is larger than that of EE, such a deformation appears near the TT singularities. 
In fact, this was the case studied in 42] which corresponds to the extreme case a — using only the TT spectrum. 
Figure O shows how the reconstructed P{k) with incorrect values of the cosmological parameters changes depending 
on a for the case that P{k) is scale invariant and that it has a peak and a dip given by Eq. (|42l) . 

In the previous case where we use only TT spectrum, even if we ended up with a spiky shape of the reconstructed 
P(fc) around some specific values of fc, it is difficult to judge if it was a real feature or simply an artifact due to 
the incorrect choice of the cosmological parameters. In contrast, in the present case where we use both TT and EE 
spectra, we can invent a new method to constrain the cosmological parameters without assuming the shape of P{k) 
from the following two observations. First, as shown in Fig. 01 in the range of 10^^ ^ ct < 10^^, P{k) is reconstructed 
rather well if we use the correct values of the cosmological parameters. On the other hand, if we use incorrect values 
in our reconstruction, we obtain a deformed shape of the reconstructed P{k). Now, as shown in Fig.|51 the location 
of the deformation in A:-space depends on the value of a we use. That is, if we take a relatively large value of a, 
say, ai ~ 10"'^^, the deformation is most prominent around the EE singularities. On the other hand, if we use a 
smaller value, say, 02 ~ 10^^, the location of the deformation moves to the TT singularities. Thus unless we adopt 
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the cosmological parameters which are close enough to the real values, the shape of Pa^ (fc), namely the reconstructed 
spectrum for a = ai, will be much different from the shape oiPa^ik)^ that for a = Thus we may obtain constraints 
on the cosmological parameters by comparison of the shapes of Pai{k) and Pa^ik). An intriguing point is that no 
assumption on the shape of P{k) is necessary to derive the constraints. 

VI. CONCLUSION 

We have derived the inversion formula including both the CMB temperature fluctuations and polarization, and 
examined the validity of our inversion method. Our method is based on a linear combination of inversion formulas 
for approximate temperature (TT) and polarization (EE) spectra in the form TT + aEE, given by Eq. H38(l . with a 
constant parameter a, combined with iteration that corrects errors in the approximate formulas. 

First we have examined the effect of the contribution of EE on the reconstructed P{k). There are singularities in 
the TT and EE inversion formulas due to the zero points of transfer functions. We have found that large numerical 
errors in the reconstructed P{k) around the singularities of the TT and EE inversion formulas, which appear when 
each formula is used independently, disappear for the choice of a '--^ 10^'*""'^^, for which the contribution of EE becomes 
comparable to that of TT. Using such an appropriate value of a, we have tested our method against several presumed 
spectral shapes of P{k). As a result, we have confirmed that the relative errors between the reconstructed P{k) and 
the original one are reduced to a few percent level after a few steps of iteration. 

We have also examined the case when the cosmological parameters are varied from the presumed real values. For the 
choice of a = 10^^, we have shown that the reconstructed P{k) is substantially deformed from its real shape near the 
EE singularities, while it is substantially deformed near the TT singularities for much smaller values of a. Since the 
errors near the singularities in the reconstructed P{k) are relatively small for any a in the range of 10^^ ^ a ^ 10^^, 
if the correct cosmological parameters are chosen, this result suggests a new method for constraining the cosmological 
parameters. Namely, by comparison of reconstructed P{k) for different values of a, we may constrain the cosmological 
parameters without any assumption on the shape of P{k). 

In this paper, we have not taken account of possible observational errors. In reality, of course, observational errors 
are by no means negligible. As shown in in the case of using the TT spectrum of the WMAP, numerical 

errors are enhanced near the singularities, which made us impossible to obtain any information of P(fc) near the 
singularities. However, our new method that substantially suppresses the numerical errors around the singularities 
may not suffer much from this difficulty. In any case, it is important to investigate the effect of observational errors on 
our new method. It is also important to investigate the possibility of using our new inversion method to constrain the 
cosmological parameters by taking account of observational errors. Since we will inevitably obtain different shapes 
of P{k) for different choices of a under the presence of observational errors, even if we knew and used the true 
cosmological parameters, it is necessary to develop a statistical criterion that tests the equivalence of the shapes of 
reconstructed P{k). It is also necessary to extend our method to include the contribution of tensor-type perturbations, 
which is pro bed by the B-mode polarization. Their contribution may be appreciable as predicted in many models of 
inflation l6ll. 163. 163. IM. l65l l6a| . These issues will be investigated in the future. 
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FIG. 1: Transfer functions f{k) {solid curve), g(k) {dashed curve), and h{k) {dash-dotted curve) for h = 0.70, Q.b ~ 0.050, SIa ~ 
0.70, Qm = 0.30 and r = 0.20. The horizontal axis is the comoving wavenumber normalized by , where d ~ 1.36 x 10'* Mpc 
is the conformal distance from the present to the LSS. f{k), g{k) and h{k) correspond to the monopole, dipole and quadrupole 
components, respectively, of the temperature anisotropy at the LSS. Since the amplitude of h{k) is much smaller than those of 
f{k) and g{k), the amplitude of h{k) is magnified by 10®. 
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FIG. 2: Comparison of the exact spectrum ''^^ {solid curves) with the approximated one '^'^'^ {dashed curves) for TT 
{top panel) and EE {bottom panel). The cosmological parameters are the same as those in Fig. Q and P{k) is assumed to be 
scale invariant. The relative errors of '^^'^ are as large as 20 — 30% for both TT and EE. 
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FIG. 3: Comparison of bf''^ between the case of a scale-invariant P{k) {solid curves), and that of P{k) with a peak and a dip 
(dashed curves), for TT {top panel) and EE [bottom panel). The cosmological parameters are the same as those in Fig.0 The 
relative differences of bf''^ between these two cases are about 5% for both TT and EE, implying the weak dependence of bf^ 
on P{k). 
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FIG. 4: Dependence of the reconstructed P{k) on the parameter a which controls the contribution of EE relative to TT, for 
the presumed form of P{k) given by Eq. (I42II in the text. The cosmological parameters are assumed to be known. The result 
of reconstruction is shown for a = 0, 10^'', 10^^, and 10^^. The horizontal axis kd roughly corresponds to £. The singularities 
f{k) = are denoted by □, and h{k) = by A. In the case of a = 10^^, for which the contributions of TT and EE are 
comparable, the errors are substantially suppressed, especially near the singularities. 
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FIG. 5: Reconstructed spectrum for P{k) with a peak and a dip, given by Eq. II42|I in the text. The top panel shows the results 
after each round of iteration, and the bottom panel shows their relative errors, \P^"^{k) — P'^'^^^ {k)\/ P'"'^^ (k) . The singularities 
f{k) — are denoted by □, and h{k) = by A. 
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FIG. 6: The same as Fig. but for P{k) with a running spectral index, given by Eq. I|48|l in the text. In this case, the errors 
are already less than 3% at the first round of iteration, and they remain stably at that level after the second and third rounds 
of iteration. 
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FIG. 7: The same as Fig. |3 but for P{k) with an osciUatory component, given by Eq. I|44|l in the text. The results are quite 
impressive because such a small oscillation like this can be reconstructed with the errors less than 3%. 
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FIG. 8: Dependence of the reconstructed P{k) on the cosmological parameters different from the presumed real values. We 
assume a scale-invariant P{k) and the cosmological parameters as ft = 0.70, ili, = 0.050, = 0.70, fim = 0.30, and r = 0.20. 
The top left panel shows the reconstruction by varying h over h = 0.60, 0.65, 0.70, 0.75, and 0.80. The top right panel shows 
the reconstruction by varying fli, over fli, = 0.040, 0.045, 0.050, 0.055, and 0.060. The bottom panel shows the reconstruction 
by varying SIa over Oa = 0.60, 0.65, 0.70, 0.75, and 0.80. In all the cases, the other parameters are fixed at the presumed 
values. Note that the value of d and the positions of the singularities depend on the cosmological parameters. Plotted are the 
singularities (/(fc) = by □ and h{k) = by A) in the case of the presumed values of the cosmological parameters. 
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FIG. 9: Dependence of the reconstructed P{k) on the parameter a in the case that one of the cosmological parameters is 
incorrect. Assumed values of the cosmological parameters are the same as those in Fig. |H] Both are the cases that we use an 
incorrect value, h — 0.75 in our reconstruction for a — 10^^ and 10^^. The top panel shows the reconstructed spectra for a 
scale-invariant P{k) and the bottom panel shows those for P{k) with a peak and a dip given by Eq. 14211 in the text. Plotted 
are the singularities in the same way as Fig. |H1 



